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Abstract 

This paper presents responses of an edge circular cantilever beam under the effect of an impact force. The beam is excited 
by a transverse triangular force impulse modulated by a harmonic motion. The Kelvin-Voigt model for the material of the 
beam is used. The cracked beam is modelled as an assembly of two sub-beams connected through a massless elastic 
rotational spring. The considered problem is investigated within the Bernoulli-Euler beam theory by using energy based 
finite element method. The system of equations of motion is derived by using Lagrange's equations. The obtained system of 
linear differential equations is reduced to a linear algebraic equation system and solved in the time domain by using 
Newmark average acceleration method. In the study, the effects of the location of crack, the depth of the crack, on the 
characteristics of the reflected waves are investigated in detail. Also, the positions of the cracks are calculated by using 
reflected waves. 
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Introduction 

Elastic wave propagation through the monitored part is of 
considerable interest in many fields. The most striking example of 
the engineering applications is detection of damage or/and 
material difference in the investigated media. By investigating 
the character of waves, the type and position of damage or/ and 
different material can be determined. 

Structural elements are subjected to destructive effects in the 
form of initial defects within the material or caused by fatigue or 
stress concentration. As a result of destructive effects, cracks occur 
in the structural elements. It is known that a crack in structure 
elements introduces a local flexibility, becomes more flexible and 
its dynamic and static behaviours will be changed. Cracks cause 
local flexibility and changes in structural stiffness. Therefore, 
understanding the mechanical behavior and the safe performance 
of edge-cracked structures are importance in designs. 

Beams can be found in many different sizes and shapes in the 
engineering applications. Circular beams are the most used and 
preferred in the machine elements because of high energy 
absorbing capability and high buckling strength. Hence, under- 
standing the mechanical behaviour circular beams are very 
important. 

In the last decades, much more attention has been given to the 
elastic wave propagation of beam structures. Teh and Huang [1] 
studied an analytical model, based on the elasticity equations, to 
investigate wave propagation in generally orthotropic beams. A 
finite element technique is developed for studying the flexural 
wave propagation in elastic Timoshenko and Bernoulli-Euler 
beams by Yokoyama and Kishida [2] . Wave propagation in a split 
beam is analyzed by treating each section separately as a 
waveguide and imposing appropriate connectivities at their joints 



by Farris and Doyle [3] . A direct mathematical approach method 
is developed to study the problem of coupled longitudinal and 
flexural wave propagation in a periodically supported infinite long 
beam by Lee and Yeen [4]. A spectral super-element model was 
used in Gopalakrishnan and Doyle [5] to model transverse crack 
in isotropic beam and the dynamic stress intensity factor was 
obtained accurately under impact type loading. Palacz and 
Krawczuk [6] investigated longitudinal wave propagation in a 
cracked rod by using the spectral element method. The use of the 
wave propagation approach combined with a genetic algorithm 
and the gradient technique for damage detection in beam-like 
structure is investigated by Krawczuk [7]. Krawczuk et al. [8] 
studied a new finite spectral element of a cracked Timoshenko 
beam for modal and elastic wave propagation analysis. Usuki and 
Maki [9] formulated an equation of motion for a beam according 
to higher-order beam theory using Reissner's principle. They 
derived the Laplace transform of the equation and investigated 
wave-propagation behavior under transverse impact. A method of 
crack detection in beam is provided by wavelet analysis of 
transient flexural wave by Tian et al. [10]. Kang et al. [1 1] applied 
the wave approach based on the reflection, transmission and 
propagation of waves to obtain the natural frequencies of finite 
curved beams. A spectral finite element with embedded transverse 
crack is developed and implemented to simulate the diagnostic 
wave scattering in composite beams with various forms of 
transverse crack by Kumar et al. [12]. The wave propagation 
model investigated herein is based on the known fact that material 
discontinuities affect the propagation of elastic waves in solids by 
Ostachowicz et al. [13]. A spectral finite element model for 
analysis of flexural- shear coupled wave propagation in laminated 
and delaminated, multilayer composite beams is presented by 
Palacz et al. [14,15]. A new spectral element is formulated to 
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analyse wave propagation in an anisotropic inhomogeneous beam 
by Ghakraborty and Gopalakrishnan [16]. Watanabe and 
Sugimoto [17] studied flexural wave propagation in a spatially 
periodic structure consisting of identical beams of finite length. 
Vinod et al. [18] investigated a formulation of an approximate 
spectral element for uniform and tapered rotating Euler-Bernoulli 
beams. Sridhar et al. [19] investigated the development of an 
effective numerical tool in the form of pseudospectral method for 
wave propagation analysis in anisotropic and inhomogeneous 
structures. An experimental method of detecting damage using the 
flexural wave propagation characteristics is proposed by Park [20] . 
Ghouvion et al. [21] studied a systematic wave propagation 
approach for the free vibration analysis of networks consisting of 
slender, straight and curved beam elements and complete rings. 
Frikha et al. [22] investigated physical analysis of the effect of axial 
load on the propagation of elastic waves in helical beams. 
Kocatiirk et al.[23] studied wave propagation of a piecewise 
homegenous cantilever beam under impact force. Kocatiirk and 
Akbas [24] investigated wave propagation of a microbeam with 
the modified couple stress theory. In a recent study, wave 
propagation and localization in periodic and randomly disordered 
periodic piezoelectric axial-bending coupled beams are studied by 
Zhu et al. [25]. Akbas, [26] studied the effect of the elastic 
foundation types on the wave propagation of the beams. 

A better understanding of the mechanism of how the crack 
effects change response of wave propagation of a circular beam is 
necessary, and is a prerequisite for further exploration and 
application of the cracked circular beams. 

In this study, wave propagation in a cantilever circular beam 
under the effect of an impact force is studied. The considered 
problem is investigated within the Bernoulli-Euler beam theory by 
using energy based finite element method. The Kelvin-Voigt 
model for the material of the beam is used. The cracked beam is 
modelled as an assembly of two sub-beams connected through a 
massless elastic rotational spring. The system of equations of 
motion is derived by using Lagrange's equations. The obtained 
system of linear differential equations is reduced to a linear 
algebraic equation system and solved in the time domain by using 
Newmark average acceleration method. The effects of the location 
of crack, the depth of the crack, on the characteristics of the 
reflected waves are investigated in detail. Also, the positions of the 
cracks are calculated by using reflected waves. 



Governing equations of intact beam 

The beam is modeled within the Euler-Bernoulli beam theory. 
According to the coordinate system (X,Y 3 %) shown in Fig. 1, based 
on Euler-Bernoulli beam theory, the axial and the transverse 
displacement field are expressed as 



u(X,Y,t)=-Y 



dv(X,t) 
dX 



v(X,Y,t) = v(X,t) 



w(X, 7,0 = 0 



(1) 



(2) 



(3) 



Where u,v and w are and z components of the displacement 
vector q, respectively, and t indicates time. 

Because the transversal surfaces of the beam is free of stress, 
then 



<5ZZ = ( 5YY = 0 



(4) 



The Kelvin-Voigt model for the material is used. The constitutive 
relations for the Kelvin-Voigt model between the stresses and 
strains become 



(5) 



where Vindicates the Young's modulus of the beam, o xx indicates 
normal stresses, £ xx indicates normal strains in the X direction, T| 
indicates the damping ratios, as follows 



(6) 



where c indicates the coefficient of damping of the beam. By using 
Eqs. (1) and (2), the strain- displacement relation can be obtained: 



du 
dX' 



,d 2 v(X,t) 
dX 2 



(7) 



Theory and Formulations 

Consider a beam of length L, diameter D, containing an edge 
crack of depth a located at a distance Li from the left end, as 
shown in Fig. 1 . One of the supports of the beam is assumed to be 
fixed and the other free. The beam is subjected to an impact force 
in the transverse direction as seen from Fig. 1. It is assumed that 
the crack is perpendicular to beam surface and always remains 
open. 



The potential energy of the beam is follows 



u,= 



EI( 



S v(X,t) 

dX 2 



)dX 



(8) 



Where / is the inertia moment of the beam. The kinetic energy of 
the beam at any instant t is 



\ * 
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\ 
\ 
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\ 



F(t) . 

\,u(x,t) 




(9) 



Where p is the mass density of the beam. The potential energy of 
the external load can be written as 



Figure 1 . A circular beam with an open edge crack subjected to 
an impact force. 

doi:1 0.1 371 /journal, pone.01 00496.g001 



F(X,t)v(X,t)dx. 



(10) 
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The dissipation function of the beam at any instant t is 



R- 



rI d 2 v(X,t)\ , 



(11) 



T=l\[pA(^]{q}) 2 ]dX 

o 



(17) 



Lagrangian functional of the problem is given as follows: 

I=T-(Ui+U e ) (12) 



R- 



[^EI(^([-^}{q})Y]dX 



(18) 



Solution method of the problem 

The considered problem is solved by using Lagrange's equations 
and time integration method of Newmark [27]. In order to apply 
the Lagrange's equations, the displacements of nodes of the 
unknown functions q (X,t) which is written for a two-node beam 
element shown in Fig. 2 are defined as follows 



{<?»}(,) = Wit) t(D vf(0 efcof 



(13) 



The displacement field of the finite element is expressed in 
terms of nodal displacements as follows 



v {e \x,t) = <p, (x)vf(t) +<p 2 W'(0 +<p 3 W'(0 +<p 4 WTO 



,(«)/ 



&(«)/• 



(14) 



"7 



where V/ and 6/ are transverse displacements and slopes at the two 
end nodes of the beam element, respectively. <Pi, <P2, 93 and 94 are 
interpolation functions and given as follows: 



The Lagrange's equations gives the following equation; 
dl d dl 



. ie) J-fTTiT) +QD k =0,k= 1,2,3, 



where 



2d. 



, k= 1,2,3,.. 



(19) 



(20) 



2^ is the generalized damping load which can be obtained from 

the dissipation function by differentiating R with respect to . 

The Lagrange's equations yield the system of equations of 
motion for the finite element and by use of usual assemblage 
procedure the following system of equations of motion for the 
whole system can be obtained as follows 



[K]{q(t)} + [D]{q(t)} + [M]{q(t)} = {F(t)} 



where 



(21) 



1^ 

[M]= | P A{<P(X)} T {<P(X)}dX (22) 



<? l (X)=l-3(X/L e ) 2 + 2(X/L e ) 3 

<? 2 (X) = L e ( - (X/L e ) + 2(X/L e f - (X/L e f) 

%(X) = 3(X/L e ) 2 -2(X/L e ) 3 

<P 4 (X) = L e ((X/L e ) 2 - (X/L e f) 



where L e is the length of the beam element. 

By substituting Equations (14) into Equations 
energy functions can be rewritten as follows: 



(15) 



(9) and (1 1), 



[Dl 



ri EI{9(X)} T {9(X)}dX 



{F(0} = 



{cp(X)} T F(X,t)dX 



x = 0 



(23) 



(24) 



,d 2 v? 



[EI{[j^]{q}Y]dX 



0, 



> 


e 











Figure 2. A two-node beam element. 

doi:1 0.1 371 /journal. pone.01 00496.g002 



(16) 



[K\- 



{9(X)} T EI{9(X)} T dX 



(25) 



where, [K] is the stiffness matrix, [D] is the damping matrix, [M] is 
mass matrix and {F(t)} is the load vector. The motion equations 
which is given by Eq. (21), are solved in the time domain by using 
Newmark average acceleration method (Newmark [27]). 

Crack modeling 

The cracked beam is modeled as an assembly of two sub-beams 
connected through a massless elastic rotational spring shown in 
Fie. 3. 
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1 -(&- 2 



Figure 3. Rotational spring model. 

doi:1 0.1 371 /journal. pone.01 00496.g003 

The bending stiffness of the cracked section kj is related to the 
flexibility G by 



1 

G 



(26) 



Cracked section's flexibility G can be derived from Broek's 
approximation (Broek [28]): 



(l-v 2 )K 2 _ M 2 dG 



2 da 



(27) 



where M is the bending moment at the cracked section, Kj is the 
stress intensity factor (SIF) under mode I bending load and is a 
function of the geometry and the loading properties as well, v 
indicates Poisson's ratio. For circular cross section, the stress 
intensity factor for Ki a single edge cracked beam specimen under 
pure bending M can be written as follow (Tada et al. [29]) 



4M h' z 



/7iaF(a/h' z ) 



(28) 



Where 



G = 



32(1 -i) 2 ) 
EttR 8 



b a z 



y(R 2 -z 2 )F 2 (a/h' z )dydz (31) 



where b and az are the boundary of the strip and the local crack 
depth respectively, are shown in Fig. 4, respectively, and written as 



b= a/R -(R-a) z 



Vr 2 -z 2 - 



(R-a) 



(32) 



(33) 



The spring connects the adjacent left and right elements and 
couples the slopes of the two beam elements at the crack location. 
In the massless spring model, the compatibility conditions enforce 
the continuities of the axial displacement, transverse deflection, 
axial force and bending moment across the crack at the cracked 
section (X = Li), that is, 



vi =v 2 , Mi =M 2 
The discontinuity in the slope is as follows: 



k T (^-^J)=k T (e 1 -e 2 )=M 1 



(34) 



(35) 



Based on the massless spring model, the stiffness matrix of the 
cracked section as follows: 



IZT, — -0.923 + 0.199(1 -sin(^f)) 4 

f <»*>= (29 » 



7ta 

cos( 2iv: ) 



Where a is crack of depth and h' z is the height of the strip, is shown 
Fig. 4, and written as 



h' z = 2VR 2 -z 2 



(30) 



where it is the radius of the cross section of the beam. 

After substituting Eq. (28) into Eq. (27) and by integrating Eq. 
(27), the flexibility coefficient of the crack section G is obtained as 



IKl 



(Cr)- 



1/G -1/G 
-1/G 1/G 



k T 
-k T 



-k T 
k T 



(36) 



The stiffness matrix of the cracked section is written according to 
the displacement vector: 



(37) 



Where 9i and 9 2 are the angles of the cracked section. With 
adding crack model, the equations of motion for the finite element 
and by use of usual assemblage procedure the following system of 
equations of motion for the whole system can be obtained as 
follows: 



A-A 

r> A r — h - — f 




Ua 

Y N/ 

Figure 4. The geometry of the cracked circular cross section. 

doi:1 0.1 371 /journal. pone.01 00496.g004 



([K] + [K] (Cr) ){q(t)} + [D]{q(t)} + [M]{q(t)} = {F(t)} (38) 



Numerical Results 

In the numerical examples, the effects of the location of crack, 
the depth of the crack, on the characteristics of the reflected waves 
are presented. In the numerical study, the physical properties of 
the beam are Young's modulus E—10 GPa, Poisson's ratio v = 0,3, 
mass density p — 2700 kg/ 'm and the damping ratio 
T| =0,0000001. The geometrical properties of the pile are length 
L=3 m and the diameter D—2 cm. The problem is analyzed 
within the framework of the Bernoulli-Euler beam theory. 
Numerical calculations in the time domain are made by using 
Newmark average acceleration method. The system of linear 
differential equations which are given by Equation (20), is reduced 
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Figure 5. The shape of the excitation impulse in the a) time domain and b) frequency domain [13]. 
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Figure 6. The relationship between first non-dimensional natural frequency and the crack depth ratio for different crack locations. 

a) Lj/L = 0.2, b) L 7 /L = 0.4. 

doi:1 0.1 371 /journal. pone.01 00496.g006 
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Figure 7. Transverse displacement at the free end of the beam, a) Intact beam, b) L 1 /L = 0A, c) L 1 /L = 03 and d) L 1 /L = 0.5. 
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Time (s) x 10 3 



Figure 8. Transverse displacement at the free end of the beam. 

D = 0.8. 

doi:1 0.1 371 /journal. pone.01 00496.g008 

to a linear algebraic system of equations by using average 
acceleration method. In the numerical calculations, the number 
of finite elements is taken as ^=100. The beam is excited by a 
transverse triangular force impulse (with a peak value 1 N) 
modulated by a harmonic function (Fig. 5) (Ostachowicz et al., 
[13]). In this study, higher frequency excitation impulse is used for 
detection of the cracks. The frequencies used in this technique are 
much higher than those typically used in modal analysis based 
methods but are lower than the frequencies used for ultrasonic 
testing. In this study, the excited frequencies lies in the range 
between 200 and 1200 kHz, with dominant one about 700 kHz. 
At such high frequencies, the response is dominated by the local 
mode and the wavelength of the excitation is small enough to 
detect incipient or potentially significant damage. (Ostachowicz 
et al., [13]). By using the equation of the wave propagation speed, 
the wave propagation speed is calculated as approximately 
0 = 4221,5 m/s. 

In order to establish the accuracy of the present formulation and 
the computer program developed by the author, the results 
obtained from the present study are compared with the available 
results in the literature. For this purpose, the first fundamental 
frequency of a cantilever circular beam with an open edge are 
calculated for different the location of crack ((Lj/L)) and the crack 
depth ratios for L = 2 m, D = 0.4 m, E = 216 GPa, p = 7850 kg/m 3 , 




i) Intact beam, b) a/D = 0.2, c) a/D = 03, d) a/D = 0.4, e) a/D = 0.6 and f) a/ 



v = 0,33 and compared with those of Kisa and Giiler [30] in the 
Figure 6. As seen from Figure 6, the present results are close to the 
results of Kisa and Giiler [30]. 

Fig. 7 illustrates the transverse displacements at the free end of 
the cantilever beam for different the crack locations (Lj/L) for the 
crack depth ratio a/D — 0.2. 

It is seen from Fig. 7 that the crack location affects significantly 
the wave propagation of the beam. It is seen from Fig. 6a that two 
waves occur (namely the excitation and the reflected wave) in the 
intact case. In Fig. 7a, the second wave occurs because of reflecting 
from the boundaries of the clamped support. In the case of the 
crack, additional secondary waves generate with first primary 
waves (see the circles) in the Fig. 7b, Fig. 7c and Fig. 7d. The 
additional secondary waves occur because of reflecting from the 
cracks. Also, it is seen from Fig. 7 that the crack locations get closer 
to the free end of the beam (namely, with the crack locations ratio 
(Lj/L) increases), the additional secondary waves appear signifi- 
cantly and the amplitude of additional secondary waves increase 
considerably. When the crack locations get closer to the clamped 
support, additional secondary waves interfere with first primary 
waves. This is because, with decrease in the distance between 
clamped support and crack, the reflected waves of the clamped 
support and the crack interfere with each other. The crack 
locations get closer to the free end, the distance between first 



PLOS ONE | www.plosone.org 



6 



June 2014 | Volume 9 | Issue 6 | el 00496 



Wave Propagation Edge Cracked Beams 




Figure 9. Transverse displacement at the free end of the beam. 

doi:1 0.1 371 /journal. pone.01 00496.g009 

primary waves and additional secondary waves increase signifi- 
cantly. 

The position of the cracks and supports can be calculated by 
using wave propagation analysis. 

For this purpose, the position of the clamped support and the 
cracks are calculated by using Fig. 7. 

Firstly, the position of the clamped support is calculated by 
using Fig. 7a: The time interval from the first wave (the excitation 
wave) and the second wave (reflected wave from the clamped 
support) is At = 7.1 x 10 _4 ,s«(see Fig. 7a). It should be noted again 
that the wave propagation speed is v— 4221.5 m/s. 

By using the relationship between velocity and time, the position 
of clamped support from the end of the beam can be calculated as 
follows; 

A r G s = vxAf = 4221.5x(7.1x 10" 4 ) = 2.9972 m 

As seen from the result, the present result is very close to the real 
position of the clamped support from the end of the beam 
{L = 3m). 

The positions of the cracks are calculated by Fig. 7b, Fig. 7c, 
and Fig. 7d for Lj/L = 0.2, Lj/L = 0.4 and Lj/L = 0.6, respectively. 

For Lj/L = 0.2. The time interval from the first wave (the 
excitation wave) and the second wave (additional wave from the 
crack) is At = 5.684 x 10 _4 .s«(see Fig. 7b). The position of the 
crack from the end of the beam is calculated as follows; 

X C rack = vxAf = 4221. 5 x (5.684 x 10" 4 ) = 2.3995m 

It is seen from the result, the present result is very close to the real 
position of the crack from the end of the beam L — 2.4 m as shown 
Fig. 7b. 

For L 1 /L = 0.4: 

The time interval from the first wave (the excitation wave) and 
the second wave (additional wave from the crack) is 




a) £=70 GPa, b) E= 100 GPa, c) £ = 300 GPa and d) £=500 GPa. 

A/ = 4.261 x 10 _4 5«(see Fig. 7c). The position of the crack from 
the end of the beam is calculated as follows; 

X C rack = vx At = 4221.5 x (4.261 x 10" 4 ) = 1.7987 m 

It is seen from the result, the present result is very close to the real 
position of the crack from the end of the beam L — 1.8 m as shown 
Fig. 7c. 

For L 1 /L = 0.6: 

The time interval from the first wave (the excitation wave) and 
the second wave (additional wave from the crack) is 
At = 2. 81925 x 10 _4 5«(see Fig. 7d). The position of the crack 
from the end of the beam is calculated as follows; 

Xcmck = v x Af = 4221.5 x (2.81925 x 10" 4 )= 1.1901 m 

It is seen from the result, the present result is very close to the real 
position of the crack from the end of the beam L — 1.2 m as shown. 

It is deduced from Fig. 7 and results that the location of the 
cracked section can be established by investigating the additional 
secondary waves. 

In Fig. 8, the transverse displacements at the free end of the 
cantilever beam for different the crack depth ratios (a/D) for the 
crack location L^L — O.b. 

As seen from Fig. 8, with the crack depth increase, the 
amplitude of additional wave increases considerably. This is 
because by increasing in the crack depth, the strength of the 
material decreases. Hence, the beam becomes more flexible. Also, 
it is seen from Fig. 8 that the generation time and location of the 
additional waves are same for different the crack depth ratios (a/ 
D). This is because; the location of the crack is not changed. It 
shows that by using wave propagation analysis, structural damages 
can be detected easily. 

In Fig. 9, the effect of the Young's modulus E on the wave 
propagation of the beam is shown for a/D = 0.3 and L^L— 0.5. 
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It is seen from Fig. 9 that Young's Modulus plays an important 
role in the additional waves. With the increased Young's Modulus 
E, the amplitude of the additional wave dramatically decreases. 
This is because by increasing in Young's Modulus, the strength of 
the material increases. Also, it is observed another result of Fig. 9 
that with the increased Young's Modulus E, the generation time 
and location of the primary and additional waves decreases. It is 
deduced from Fig. 9 that Young's Modulus is very effective for 
reducing the negative influence of the cracks. 

Conclusions 

Wave propagation in an edge circular cantilever beam under 
the effect of an impact force is investigated. The effects of the 
location of crack, the depth of the crack and Young's Modulus on 
the wave propagations of the circular beam are investigated in 
detail. The following conclusions are reached from the obtained 
results: 

References 

1. Teh KK, Huang GG (1981) Wave propagation in generally orthotropic beams. 
Fibre Science and Technology 14: 301 310. 

2. Yokoyama T, Kishida K (1982) Finite element analysis of flexural wave 
propagation in elastic beams. Technology Reports of the Osaka University 32: 
103-112. 

3. Farris TN, Doyle JF (1989) Wave propagation in a split Timoshenko beam. 
Journal of Sound and Vibration 130: 137 147. 

4. Lee SY, Yeen WF (1990) Free coupled longitudinal and flexural waves of a 
periodically supported beam. Journal of Sound and Vibration 142: 203 211. 

5. Gopalakrishnan S, Doyle JF (1995) Spectral super-elements for wave 
propagation in structures with local non-uniformities. Gomput. Methods Appl 
Mech Eng 121:77-90. 

6. Palacz M, Krawczuk M (2002) Analysis of longitudinal wave propagation in a 
cracked rod by the spectral element method. Computers and Structures 80: 
1809-1816. 

7. Krawczuk M (2002) Application of spectral beam finite element with a crack and 
iterative search technique to damage detection. Finite Element In Analysis and 
Design 38(6):537-548. 

8. Krawczuk M, Palacz M, Ostachowicz W (2002) The Dynamic Analysis of a 
Cracked Timoshenko Beam By The Spectral Element Method. Journal of 
Sound and Vibration 264: 1139-1153. 

9. Usuki T, Maki A (2003) Behavior of beams under transverse impact according to 
higher-order beam theory. International Journal of Solids and Structures 40: 
3737-3785. 

10. Tian J, Li Z, Su X (2003) Crack detection in beams by wavelet analysis of 
transient flexural waves. Journal of Sound and Vibration 261: 715—727. 

1 1 . Kang B, Riedel CH, Tan CA (2003) Free vibration analysis of planar curved 
beams by wave propagation. Journal of Sound and Vibration 260: 19—44. 

12. Kumar DS, Mahapatra DR, Gopalakrishnan S (2004) A spectral finite element 
for Wave propagation and structural diagnostic analysis of composite beam with 
transverse crack. Finite Elements in Analysis and Design 40: 1729-1751. 

13. Ostachowicz W, Krawczuk M, Gartmell M, Gilchrist M (2004) Wave 
propagation in delaminated beam. Computers and Structures 82: 475 483. 

14. Palacz M, Krawczuk M, Ostachowicz W (2005) The Spectral finite element 
model for analysis of flexural-shear coupled wave propagation, part 1: laminated 
multilayer composite. Composite Structures 68: 37-44. 

15. Palacz M, Krawczuk M, Ostachowicz W (2005) The spectral finite element 
model for analysis of flexural-shear coupled wave propagation. Part 2: 
Delaminated multilayer composite beam. Composite Structures 68: 45-51. 



(1) The wave propagation analysis can easily be used for crack 
detection procedures within structures. 

(2) The crack locations can be established by investigating in the 
additional waves. 

(3) The magnitude of the crack is very effective in the amplitude 
of additional waves. 

(4) The crack locations and the crack depth have a great 
influence on the wave propagation of the circular beam. 

(5) There are significant differences of the wave propagation for 
the cracked and intact circular beams. 

(6) Young's Modulus is very effective for reducing the negative 
influence of the cracks. 



Author Contributions 

Conceived and designed the experiments: SDA. Performed the experi- 
ments: SDA. Analyzed the data: SDA. Contributed reagents/materials/ 
analysis tools: SDA. Wrote the paper: SDA. 



16. Chakraborty A, Gopalakrishnan S (2005) A spectral finite element for axial- 
flexural-shear coupled wave propagation analysis in lengthwise graded beam. 
Computational Mechanics 36: 1-12. 

17. Watanabe Y, Sugimoto N (2005) Flexural wave propagation in a spatially 
periodic structure of articulated beams. Wave Motion 42: 155 167. 

18. Vinod KG, Gopalakrishnan S, Ganguli R (2007) Free vibration and wave 
propagation analysis of uniform and tapered rotating beams using spectrally 
formulated finite element. International Journal of Solids and Structures 44: 
5875-5893. 

19. Sridhar R, Chakraborty A, Gopalakrishnan S (2007) Wave propagation analysis 
in anisotropic and inhomogeneous uncracked and cracked structures using 
pseudospectral finite element method. International Journal of Solids and 
Structures 43: 4997-5031. 

20. Park J (2008) Identification of damage in beam structures using flexural wave 
propagation characteristics. Jounal of Sound and Vibration 318: 820 829. 

21. Chouvion B, Fox CHJ, Mc William S, Popov AA (2010) In-plane free vibration 
analysis of combined ring-beam structural systems by wave propagation. Jounal 
of Sound and Vibration 329: 5087-5104. 

22. Frikha A, Treyssede F, Cartraud P (201 1) Effect of axial load on the propagation 
of elastic waves in helical beams. Wave Motion 48: 83—92. 

23. Kocatiirk T, Eskin A, Akba§ §D (2011) Wave propagation in a piecewise 
homegenous cantilever beam under impact force. International Journal of the 
Physical Sciences 6: 4013-4020. 

24. Kocatiirk T, Akbas, §D (2013) Wave propagation in a microbeam based on the 
modified couple stress theory. Structural Engineering and Mechanics 46: 417- 
431. 

25. Zhu H, Ding L, Yin T (2013) Wave propagation and localization in a randomly 
disordered periodic piezoelectric axial-bending coupled beam. Advances in 
Structural Engineering 16: 1513-1522. 

26. Akbas, §D (2014) Wave propagation analysis of edge cracked beams resting on 
elastic foundation. International Journal of Engineering & Applied Sciences 6: 
40-52. 

27. Newmark (1959) A method of computation for structural dynamics. ASCE 
Engineering Mechanics Division 85: 67—94. 

28. Broek D (1986) Elementary engineering fracture mechanics. Dordrecht: 
Martinus Nijhoff Publishers. 

29. Tada H, Paris PC, Irwin GR (1985) The Stress Analysis of Cracks Handbook. 
St.Louis: Paris Productions & (Del Research Corp.). 

30. Kisa A, Giiler MA (2007) Free vibration analysis of uniform and stepped cracked 
beams with circular cross sections. International Journal of Engineering Science 
45: 364-380. 



PLOS ONE | www.plosone.org 



8 



June 2014 | Volume 9 | Issue 6 | el 00496 



